Numerical Simulation of Three-Dimensional Dendrites using Coupled Map Lattices 
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Three dimensional dendrites are studied with a coupled map lattice model. We study the fractal 
, dimensions, the /(a) spectrum, the size distribution of sidebranches, and the envelope formed by 

sidebranches. 
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Branched patterns such as diffusion-limited aggregation (DLA) and dendrite appear in various diffusion fields under 
[ strongly nonequilibrium conditions Dendrites with rough surfaces have been intensively studied by theoretical 

' I analyses, numerical simulations, and experiments 0,0. Two-dimensional DLA and dendrites were studied fairly well 
,—1 theoretically and experimentally; however, three-dimensional dendrites have not been sufhciently studied. Huang and 
, Glicksman studied dendritic tips and sidebranch structures of succinonitrilej6i, i7|]. They found the law vp^ =const for 
^ ' a dendritic tip, where v is the tip velocity and p is the radius of curvature, and found that the spacing of neighboring 

a dominant sidebranches increases as Az ~ z^'^, where z is the distance from the dendritic tip and Az is the spacing 
between neighboring dominant sidebranches. The dendritic tip is well approximated by a paraboloid and the cross 
-4-^ ^ section is almost circular (axisymmetric) for succinonitrile. More detailed studies were performed in microgravity 
^ • using images taken from the space shuttle [8| . Sidebranches can grow owing to thermal noises [§| and develop far from 
I the dendritic tip. Li and Beckermann found various scaling laws for sidebranches in an experiment on succinonitrile 
(lOj . The envelope formed by sidebranches obeys X oc z'^-^^g ^^^^ ^j^g volume V between the dendritic tip and z obeys 
V oc z^'^. The cross section of the dendritic tip is not always axisymmetric and a fin structure develops in some 
dendrites. Brener and Temkin studied such non-axisymmetric three-dimensional dendrites theoretically and found 
that the shape of the dendritic tip can be described by z = a|x|^/^, where x is the length scale of the fin structure 
rH ' perpendicular to the growth direction. It is different from the parabolic shape z = a\x P Bisang and Bil gram 

Q ' confirmed the scaling law of z = a|a;|^/'^ in an experiment on xenon dendrites [isj . Wittwer and Bilgram measured 
O ' the envelope of sidebranches of xenon dendrites and estimated X ~ 0.5 • z [T^ . 

On the other hand, we proposed a coupled map lattice model as a simple simulation method of generating various 
' growth patterns such as DLA, DBM, and dendrites [H, [l^. The envelope shape of two-dimensional dendrites was 
^ , studied using a coupled map lattice ^3,- The height distribution of sidebranches was studied using a coupled map 
T-H lattice and compared with a needle model [l^- The height distribution obeys the power law p{h) ^ h~^-^. The 
power-law distribution of sidebranches was confirmed in an experiment on quasi- two-dimensional NH4CI crystals [l9| . 
■ In this study, we construct three-dimensional dendrites using a coupled map lattice, and study the fractal dimension, 
/(a) spectrum, envelope shape, and height distribution of the sidebranhches of three-dimensional dendrites. 
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^ ' II. THREE-DIMENSIONAL COUPLED MAP LATTICE 

^ ^ , Dendritic crystals grow from a supersaturated solution. Assuming the cubic symmetry for such crystals, a coupled 

k> map lattice model on a cubic lattice is proposed. The coupled map lattice for the solution growth is composed of two 

5_j processes. One is the diffusion of the solution: 



u'^{i,j, k) = Un{i,j, k) + D{un{i + l,j,k) + Un{i - l,j,k) + w„(i,j + 1, k) 

+Un{i,j -l,k) + Un{i,j,k + 1) + Un{iJ, fc - 1) - 6u„(z, j, fc)}, (1) 

where u, and u' denote the dimensionless concentration of the solution, n is the number of steps, (i, j, k) denotes the 
lattice point, and D is the diffusion constant. The second step is the growth at the interface. The order parameter 
m(i, j, k) is introduced at each lattice point to express the growth. The order parameter m(i,j, k) is set to be at sites 
of the solution, and 1 at sites of the crystals. The order parameter m{i,j,k) changes only at the interface between 
solution sites and crystal sites. The growth rules of m and u at interface sites are written as 

m„+i(i,j, fc) = mn(i,j,k) + c{u[^{i,j,k)-Uc), 

Un+i{i,j,k) = uji,j,k) - c{u'„{i,j,k) -Uc), (2) 
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FIG. 1: (a) 3D plot of a three-dimensional dendrite for ito = 0.01. A region for fc > 60 is plotted, (b) Cross section of the 
three-dimensional dendrite at fc = 0. 

(a) Cb) 

I 1 100000 E . . — —5 




100 200 300 10 100 

i R 



FIG. 2: (a) Cross section of a three-dimensional dendrite for u = 0.01 and A = 0.04. (b) Relation of the radius R of gyration 
and the number of crystal sites, S. The dashed line denotes S ~ E?'^^ . 

where Uc denotes the equilibrium concentration of the solution, and c denotes the kinetic coefficient of deposition 
from the solution into the crystals at interface sites, which is proportional to the local supersaturation u — Uc- The 
conservation law of mass is satisfied in our coupled map lattice model during growth at the interface, that is, the 
concentration of the solution around crystal sites decreases as crystalization proceeds. If m„+i(i,j, fc) goes over a 
threshold nic ^ 1, the (i, j, k) site is set to be a site of crystals. The same processes are repeated. The threshold rric is 
1 for the simplest model. Randomness can be incorporated by setting the threshold rric as a random number between 
1 — A < rric < 1 + A. We assume that Uc is when the effect of surface tension is not considered. The effect of surface 
tension is incorporated in our model by counting the number of neighboring sites of crystals, N , for each interface 
site, and setting Uc — ^{N — 9). Here, the neighboring sites are sites denoted by (i ± 1, j ± 1, fc ± 1) around the site 
(j, J, fc). For a flat surface, iV = 9 and Uc = 0. For a tip site, TV = 1 and Uc = —87 < 0; then growth is suppressed at 
such a protruded site. 

III. FRACTAL DIMENSION AND /(a) SPECTRUM 

We have performed numerical simulations of a coupled map lattice. A typical three-dimensional dendrite is shown 
in Fig. 1. A cubic lattice of [—L/2, L/2] x [— L/2, L/2] x [— L/2, L/2] is used in numerical simulation. Fixed boundary 
conditions are imposed as u{i, j, k) — uq at i — ±L/2, j = ±L/2 or fc = ±L/2. A point seed is set at the center (0, 0, 0) 
where L = 300. The initial concentration is u{i,j,k) = uq — 0.01. The other parameters are c = 1,7 = 0,D = 0.15, 
and A = 0. Figure 1(a) shows a 3D plot of the dendrite for fc > 60, and Fig. 1(b) shows the cross section at fc = 0. 

Figure 2(a) shows the cross section of a larger three-dimensional dendrite for i > and j > at fc = . The 
parameter values are uq = 0.01,7 = 0,L — 700, and A = 0.04. During the dendritic growth, the total number of 
crystal sites, S'„, is counted at a time step n, and the radius of gyration, Rn, at the same time step n is calculated 
using i?^ = J2i j k '"'J j k/'^n, where r^j-^fc is the distance between a crystal site fc) and the origin. Figure 2(b) shows 
the relation of i?„ and Sn- The power law relation S ^ with D ~ 2.25 is observed. The exponent D was evaluated 
using the data for R > 30. D is the fractal dimension of a three-dimensional dendrite. The fractal dimension of the 
three-dimensional DLA is D ^ 2.5, but there are few reports on the fractal dimension of three-dimensional dendrites. 
On the other hand, the fractal dimension of two-dimensional dendrites is D ^ 1.5, which has also been checked using 
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FIG. 3: (a) Relation of q and Dg. (b) Relation of a and /(q) for u = 0.01, 7 = 0, and A = 0.04. 

our two-dimensional coupled map lattice. 

Next we calculated the f{a) spectrum for the growth probability [20^. We considered the fact that the concentration 
u{i,j,k) at the interface site is proportional to growth probability. The partition function Zq(L) is calculated as 



Z,(L) = 5]{P(z,j,fc)}^ 



(3) 



where P{i, j, k) = u{i,j, fc)/ ^ u{i, j, fc), and the summation is taken for all the interface sites. The size of the dendrite, 
L„, grows with time; here L„ is measured as the distance between the tip position and the origin. The generalized 
dimension Dq is calculated using the partition function as 



-1 



q - 1 L„ 



lim 



\nZq{Lr. 
lnL„ 



(4) 



Numerically, limL^^oojln Zq(L„)}/ lni„ was evaluated at the linear slope of the plot of In Zq{Ln) vs lnX„ for X„ > 30. 
The f{a) spectrum is calculated using 



= ^{{q-l)Dq}, 
fiaiq)) = qa{q)-{q-l)Dq 



(5) 



Figure 3 shows the Dq spectrum and f{a) spectrum of the dendrite at uq = 0.01,7 = 0,i = 700, and A = 0.04. 
The peak value of /(a) is 2.24, which is the fractal dimension Dq. The information dimension Di is evaluated as 
Di — 1.87, and I? 00 is evaluated as D^o — 1.34. The minimum a is evaluated as amin — 1-24. The f{a) spectrum is 
rather wide in comparison with that of the two-dimensional dendrite. The spectrum for D{q) for q < —2.4 could not 
be accurately calculated, because a clear linear relation of lnZ^(L„) and lnL„ was not obtained or the data scattered 
widely for small values of q. Then, the /(a) spectrum is not accurately calculated for a > 56. We have also calculated 
the /(a) spectrum of dendrites at other parameter values such as uq = 0.005,0.02 and A = 0,0.02, but the /(a) 
spectrum does not change markedly. 



IV. ENVELOPE OF SIDEBRANCHES 



We have studied three-dimensional dendrites with surface tension using the coupled map lattice model by incor- 
porating the parameter 7. Figure 4(a) shows a cross section of a dendrite at fc = and i > 50 for c = 0.8, 
7 = 4 X 10~^, Mo = 0.005, and A = 0. A smooth interface appears owing to the surface tension effect. The interface is 
approximated as j = 27 — 0.064i for i < 200. Figure 4(b) shows the magnification of the tip region. Near the tip region, 
the interface is approximated by j = 0.77 • (328 — i)°-^. This is consistent with the scaling law x ^ by Brener, 
because z is the distance from the dendritic tip and x is the height of the fin structure in the notation by Brener, and 
z and X respectively correspond to 328 — i and j in our coupled map lattice. Figure 4(c) shows a cross section in the 
(j, k) space at i = 50. A two-dimensional dendrite appears in this cross section. The ratio of the tip velocity in the 
j-direction of this two-dimensional dendrite and the tip velocity in the i-direction of the three-dimensional dendrite 
shown in Fig. 4(a) is equal to the slope 0.064 of the dashed line in Fig. 4(a). 

The smooth interface becomes unstable for 7 < 2.7 x 10~^. Figure 5 shows snapshot patterns for 7 — 2.5 x 
10~^, 1.5 X 10~^, and for uq = 0.005 and A = 0. Sidebranches appear in a sporadic manner far from the dendritic 
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FIG. 4: (a) Cross section of a three-dimensional dendrite in the space at A; = for i > 50. The parameter 7 is 4 x 10~^. 
The dashed line is j — 27 — 0.064i. (b) Cross section of a three-dimensional dendrite in the {i,j) space for i > 250. The dashed 
curve is j = 0.77 • (328 — j)" ''. (c) Cross section of a three-dimensional dendrite in the (j, k) space at i = 50. 
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FIG. 5: Cross section of a three-dimensional dendrite in the space at fc = for i > 50. The parameter 7 is (a) 2.5 x 10 
(b) 1.5 X 10-^ and (c) 0. 



tip at 7 = 2.5 X 10^^. The number density of sidebranches increases as 7 decreases. We have constructed an envelope 
of sidebranches along the i axis by connecting the tip positions jt{i) of sidebranches that are taller than any tip 
positions jt(j') at i', satisfying i' > i. We have performed many numerical simulations and constructed an average 
shape of the envelope of sidebranches. In these numerical simulations, the random threshold of A = 0.04 is assumed 
only on the three lines of (i, 0, 0), (0, j, 0), and (0,0, fc) passing through the origin. Figure 6(a) shows a snapshot 
pattern for uq = 0.005 and 7 = 0. The dashed Une is j = 0.39 • (440 — i). Figure 6(b) shows average shapes of the 
envelopes of sidebranches at 7 = and 7 = 1.5 x 10^^, where z is the distance from the dendritic tip at j = fc = 0, 
and X is the height of the envelope at z. The envelopes are approximated at X ~ 0.39 • z for 7 = and X ~ 0.37 ■ z 
for 7 = 1.5 X 10~^ for z < 150, and the envelopes bend slightly downward as X ~ 2;°'^ for z > 150. These results are 
consistent with the experimental results. 




FIG. 6: (a) Snapshot pattern at uo = 0.005 and 7 = 0. (b) Envelopes of sidebranches for 7 = and 7 = 1.5 x 10 
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FIG. 7: (a) Competition of two branches at ug — 0.005 The system size is 22 x 280 x 280. (b) Time evolution of the two tip 
positions at i = 7 and 18. (c) Relation of vi and the difference between the two tip positions for I — 11 and 19. 
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FIG. 8: (a) Snapshot pattern at uo = 0.005 and 7 = grown from a linear seed at j — 0. (b) Size distribution p{h) of branch 
size h. The dashed line denotes a power law of exponent 1.35. 

V. COMPETITIVE GROWTH AMONG SIDEBRANCHES 

Sidebranches compete with each other via the diffusion field. When a sidebranch grows slightly faster than the 
neighboring sidebranches, the taller sidebranch grows even faster because the sidebranch is in contact with a denser 
diffusion field. In a previous paper, we performed a numerical simulation of the competitive growth of two sidebranches 
using a two-dimensional coupled map lattice ^1 7,]. We have found that the growth velocity of a shorter sidebranch 
decays in an exponential manner, and that decay rate decreases with the power law of the distance between the two 
sidebranches, I. We have performed a similar numerical simulation in three dimensions. We have used a rectangular 
box oi Lx y< L X L for numerical simulation, and periodic boundary conditions are imposed in the cc-direction. As 
an initial condition, we have set two seeds of crystals of different sizes at i = ii — Lx/2 — (/ + l)/2 and i = 12 = 
Lx/2 + + l)/2. Figure 7(a) shows a snapshot at n = 150000 for uq = 0.005,7 = 0, A = 0,L = 280, = 22, and 
I — 11. The right branch at i = 17 grows faster and wins the competition, because the initial seed is slightly larger. 
Figure 7(b) shows the time evolution of the tip positions jti and jt2- The tip position jti of the left branch stops to 
grow because the right branch dominates the diffusion field. Figure 7(c) shows a semilogarithmic plot of the growth 
velocity vi of the left branch as a function of the difference A — jt2 ~ jti of the tip positions for Z = 11 and 19. The 
velocity vi decays slower than an exponential decay. The dashed curves are fitting curves: 0.0012 exp(—0.75A"-^^) for 
I = 11 and 0.0004 exp(—0.15A°-^^) for / = 19. It is different from the exponential decay observed in two-dimensional 
dendrites in our previous paper. We do not understand the reason for the stretched exponential decay, however, the 
slower decay implies that the suppression effect by the longer branch is weaker in three dimensions. This is probably 
because the shorter branch can grow slightly owing to the diffusion of u in the fc-direction. 

Next, we have performed a numerical simulation from a linear seed at j = and fc = to study the competitive 
growth of many sidebranches. Figure 8(a) shows a snapshot pattern at fc = for uq = 0.005,7 = 0: and A = 0. 
The initial m{i,j, k) at interface sites around the linear seed is a random number between and 0.2. Many branches 
grow and complete with each other. Several large branches survive and grow. Finer-scale sidebranches develop on 
the surviving large branches. The growth of other small branches is suppressed owing to the larger sidebranches. 
Figure 8(b) shows the size distribution p{h) of branches on a double-logarithmic scale. The size distribution was 
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constructed from 22 branch patterns in the i — j space and i — k space obtained from 11 random initial conditions. 
The approximate curve is the power law of p(ft-) ~ l/Zi^ "^^. The exponent of this power law was 1.9 in two-dimensional 
coupled map lattices. The exponent in our three-dimensional coupled map lattice is rather small, which might be 
related to the fact in Fig. 7 that the competition is weak in three dimensions. 



VI. CONCLUSIONS 



We have performed numerical simulations of three-dimensional dendrites using coupled map lattices. We have 
succeeded in generating well-developed three-dimensional dendrites. We have calculated the fractal dimension and 
/(a) spectrum. We have checked the scaling law x ~ z^^^ of the smooth tip region by Brener by incorporating a 
surface tension effect, and found a transition from a smooth surface to a well-developed dendrite. Furthermore, we 
have studied competitive dynamics among sidebranches, and found some difference of what from two-dimensional 
dendrites. 
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